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1.  Introduction 


Understanding  the  scattering  behaviors  of  trees  plays  a  critical  role  in  determining 
the  performance  of  foliage-penetration  radar  systems.  The  foliage-penetration  radar 
modality  of  interest  in  this  study  consists  of  an  airborne  sensing  system  that 
illuminates  a  scene  with  a  low  frequency  signal  (typically,  at  low-grazing  angles), 
and  then  measures  the  backscattered  responses  in  the  far  field.  Subsequently,  the 
detection  and  surveillance  of  ground  targets  obscured  by  the  trees  can  be  carried 
out  with  various  clutter  cancellation,  moving  target  indication,  and  synthetic 
aperture  radar  imaging  processing  methods.  In  the  low  frequency  region,  in  terms 
of  the  scattering  effects,  the  contributions  from  the  leaves  and  small  branches  are 
often  assumed  to  be  negligible;  and  while  the  inclusion  of  the  contributions  from 
large  branches  can  modify  the  overall  response  to  an  extent,  the  most  salient 
scattering  features  in  the  backscattering  direction  are  generated  primarily  by  the 
tree  trunks.'  Characterizing  the  electromagnetic  responses  of  the  tree  trunks,  then, 
constitutes  an  integral  part  of  the  detection  algorithm  development  and  radar 
performance  evaluation  process. 

Although  full-wave  simulations  of  large  forest  scenes  have  been  carried  out  in 
previous  studies,'’^  numerical  solutions  are  rather  computationally  intensive.  An 
analytical  approach  is  not  only  advantageous  in  computational  efficiency,  but  also 
in  providing  a  more  satisfying  physical  interpretation  of  the  scattering  effects.  In 
accordance  to  the  discrete  scatterer  approach,^  a  tree  trunk  can  be  treated  as  a  finite- 
length,  dielectric  circular  cylinder,  the  scattering  solution  of  which  can  be  obtained 
with  approximate  analytical  methods  such  as  surface  current  integration  or 
volumetric  current  integration."'^^  (Note  that  a  physical  optics  technique  can  also  be 
applied,  but  such  an  approach  is  more  appropriate  at  the  higher  frequency  region. 
For  the  range  of  parameters  of  interest  here,  the  radii  of  curvature  of  the  trunks  are 
not  always  large  compared  to  the  wavelength.)  The  surface  current  integration 
method  exploits  the  supposition  that  the  surface  fields  of  the  finite-length  cylinder 
are  the  same  as  those  of  the  infinite-length  case,  whereas  the  volumetric  current 
integration  method  makes  a  similar  assumption  but  for  the  internal  fields."'’^  With 
the  stated  approximation,  the  scattering  matrix  elements  of  the  finite-length 
structure  can  then  be  calculated  using  the  equivalence  principle.  The  validity  of 
both  methods  has  been  investigated  in  previous  studies  for  scatterers  in  free 

7  9  10 

Space.  ’  ’ 

The  focus  of  the  current  work  is  on  the  volumetric  current-based  approach,  since  it 
has  been  shown  that  this  method  is  more  accurate  relative  to  the  surface  current- 
based  one.’  Detailed  derivations  of  the  relevant  formulations  are  discussed  in  Part 
I  of  this  study.^  However,  a  comprehensive  evaluation  of  the  accuracy  of  this 
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method  for  scatterers  located  above  a  dielectric  ground  has  not  been  presented.  As 
such,  this  work  is  organized  as  follows.  In  Section  2,  an  overview  of  the  analytical 
formulations  for  the  scattering  of  a  tree  trunk  located  above  a  dielectric  ground  is 
put  forth.  In  Section  3,  the  region  of  validity  of  the  closed-form  backscattering 
solution  as  functions  of  trunk  length,  trunk  radius,  incidence  angle,  polarization, 
and  trunk  tapering  is  deduced  by  comparing  the  results  with  those  from  a  full-wave 
solution;  the  responses  are  also  examined  in  the  imaging  domain  to  elucidate  the 
relevancy  of  various  scattering  mechanisms.  Finally,  in  Section  4,  a  summary  of 
the  work  is  given. 


2.  Analytical  Solution 

As  discussed  in  Part  I  of  this  study, ^  a  particular  closed-form  solution  for  the 
scattering  from  the  finite-length,  dielectric  circular  cylinder  (with  length  L  and 
radius  a)  in  free  space  can  be  obtained  by  approximating  the  internal  fields  of  the 
cylinder  as  those  of  the  infinite-length  case.  A  scattering  matrix  can  be  subsequently 
constructed  with  an  application  of  the  volumetric  equivalence  principle  by 
propagating  the  fields  radiated  by  the  electric  polarization  current  inside  the 
cylindrical  structure.  As  pointed  out  by  Lang  and  de  Matthaeis,’  '**  such  a  solution 
is  reasonably  accurate  in  the  forward  scattering  cone,  and  the  accuracy  improves 
with  increasing  cylinder  length,  for  a  given  radius.  Note  that  by  integrating  over  a 
finite  cylindrical  volume,  this  method  somewhat  implicitly  includes  the  scattering 
from  the  ends  of  the  finite-length  structure;  however,  this  scattering  component 
may  not  be  captured  accurately.  Here,  the  overall  solution  above  is  referred  to  as 
the  volumetric  current  integration  method,  and  the  relevant  formulations  from  Part 
I  of  this  study  are  tabulated  in  the  Appendix.  For  the  case  when  the  tree  trunk  is 
located  above  a  finite-conducting  ground  (with  the  air-ground  interface  at  z  =  0),  a 
total  scattering  matrix  can  be  constructed  with  a  multi-ray  approach  by  coherently 
summing  the  direct  wave  and  the  ground-reflected  waves.  Succinctly,  the  total 
scattering  matrix — expanded  into  its  first  4  main  contributions — can  be  written 


(4,, i, '■kr(t4./(4^.4)T(i,)e 
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where  ^  is  the  location  of  the  center  of  the  scatterer;  -2(r^ -7)2 ; 

k.^  =k.-2 {k^  ■z'jz;  =  4  - 2 (4  •  -z) ^  (•’ •)  cylinder  free  space  response 

given  in  the  Appendix;  and  r  ^ )  is  the  ground  reflection  coefficient  matrix 


r 
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with  ^  j  and  ^  j  as  the  ordinary  horizontally  and  vertically  polarized 

ground  reflection  coefficients,  respectively.  The  incident  and  scattered  wave 
directions  ( k^  and  ,  respectively)  and  their  associated  polarization  vectors  are  as 

defined  in  the  Appendix.  The  first  term  in  Eq.  1  is  the  direct  wave  from  the  cylinder; 
the  second  and  third  terms  correspond  to  the  single  ground-bounce  returns;  and  the 
last  term  is  the  double  ground-bounce  contribution.  Note  that  in  the  backscattering 

direction  {k^- -k. ),  the  dominant  response  comes  from  the  single  ground-bounce 

returns — ^that  is,  from  the  dihedral  effect  of  the  trunk  and  the  ground  surface.  The 
formulation  in  Eq.  1  essentially  exploits  3  approximations:  first,  the  trunk 
volumetric  current  is  estimated  as  that  from  the  infinite-length  cylinder  in  free 
space;  second,  the  ground  effects  on  that  current  are  taken  as  negligible;  and  finally, 
the  scattering  in  the  presence  of  the  ground  is  heuristically  captured  as  rays. 

In  terms  of  modeling  the  physical  geometry,  the  above  developments  assume  the 
tree  trunk  is  a  circular  cylindrical  structure  with  a  constant  radius.  A  more  realistic 
representation  is  a  cylindrical  structure  that  has  a  varying  radius  along  its  length.  In 
fact,  most  tree  trunks  in  nature  have  a  tapered  profile  defined  by  an  exponential 
expansion  (or  flare)  near  the  base."  Accordingly,  circular  cylinders  with  and 
without  tapering  along  their  lengths  are  considered  in  this  work  (Fig.  1).  Scattering 
from  a  tapered  trunk  is  treated  by  first  segmenting  the  structure  into  discrete 
constant-radius  subsections  along  its  length;  then,  the  response  of  each  subsection 
is  deduced  using  Eq.  1 ;  and  finally,  the  total  response  of  the  whole  structure  is 
calculated  by  summing  the  responses  of  all  the  subsections. 
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Fig.  1  Truuk  proOles^  considered  in  this  study:  a)  uu tapered  cii'cular  cyliudiical  trunk, 
b)  linearly  tapered  circular  cylindrical  trunk,  and  c)  nonlinear ly  tapered  circular  cylindiical 
trunk*  Each  structure  has  base-to-apex  height  (or  length)  L  and  base  radius  a. 


3.  Validation  with  Full-Wave  Solution 


The  accmacy  of  the  analytical  scattering  solution  for  a  tree  trunk  located  above  a 
ground  is  investigated  by  comparing  the  backscattering  results  with  those  from  a 
finite-difference  time-domain  (FDTD)  solver.  Only  the  co-polarized  retiuns 
are  considered  in  this  work.  (Formally,  the  backscattered  cross-polarized  responses 
are  zero  for  the  stnictmes  studied  here.)  hi  the  frequency  domain,  the  average  error 
in  the  closed-fonn  solution  is  examined  as  functions  of  the  tmnk  length,  trunk 
radius,  and  incidence  angle  for  3  different  stmctures — ^the  untapered  circular 
cylindrical  hunk,  the  linearly  tapered  circular  cylindrical  hunk,  and  the  nonlineaiiy 
tapered  circular  cylindiical  hunk.  The  average  eiror  is  defined  as 


where  5^^  denotes  the  reference  (FDTD)  solution,  and  the  averaging  operation  is 
perfoimed  over  the  frequency  band  Af  =  [200  MHz,  500  MHz]  (P-band)  with 
=  301  frequency  points.  Trunk  profiles  with  dimensions  hi  the  range  of 
o.  lA  <  i  <  lOA  and  o.  lA  <  <7  <  A  are  treated  in  steps  of  0.  ^  where  Xc  is  the 

free  space  wavelength  at  the  center  of  the  band,  and  are  the  normalized 

length  and  radius  of  the  structure,  respectively.  That  is,  1000  different  geometries 
are  simulated  altogether — for  each  of  the  3  tapering  profiles.  Note  that  some 
cases — for  example,  those  with  small  length-to-radius  ratios — ^aie  defined  by  a 


combination  of  parameters  that  is  unlikely  to  occur  in  nature  and  therefore  are  not 
realistic;  nevertheless,  these  cases  are  included  in  the  analysis  for  the  sake  of 
completeness.  In  this  study,  the  trunk  and  the  ground  assume  the  following 
dielectric  properties:  the  trunk  has  a  homogeneous  dielectric  composition  with  a 
real  part  for  the  relative  dielectric  constant  of  ^  =13.90  and  a  conductivity  of 

at  =  39  mS/m  (specifically,  these  values  correspond  to  an  aspen  tree  trunk  with  55% 
moisture  content,  as  indicated  by  the  measurement  data  appearing  in  Koubaa  et 
al.);'^  and  the  ground  has  g  =  5.45  and  ag  =  20  mS/m  (these  values  correspond 
to  a  soil  with  8%  clay,  30%  sand,  62%  silt,  and  10%  moisture  content). 

3.1  Untapered  Circular  Cylindrical  Trunk 

The  first  case  analyzed  is  the  regular  circular  cylinder  trunk  model  (Fig.  la). 
Figures  2  and  3  show  the  average  error  in  the  analytical  solution  as  referenced  to 
the  FDTD  solution  as  functions  of  the  trunk  length,  trunk  radius,  and  incidence 
angle  for  the  vv  and  hh  responses.  In  order  to  limit  the  dynamic  range  of  these  error 
maps,  error  data  points  greater  than  15  dB  are  simply  set  to  15  dB.  A  constant  error 
contour  is  also  extracted  from  each  error  map  and  is  plotted  in  Figs.  4  and  5;  here, 
the  3-dB  error  line  is  used,  that  is,  if  the  trunk  length  and  radius  values  lie  in  the 
region  above  this  line,  then  the  average  error  is  expected  to  be  less  than  3  dB.  (Note 
that  for  the  vv  response  at  Ot  =  15°,  as  the  error  does  not  fall  below  3  dB  for  the 
region  >  0.24  over  the  range  of  length  values  simulated,  only  a  partial  error 

line  is  drawn.  For  this  special  case,  additional  analysis  of  trunks  with  lengths  up  to 
20Xc  is  carried  out;  however,  a  complete  error  line  still  cannot  be  determined.  The 
reasons  for  the  poor  accuracy  are  explained  below.) 
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Fig.  2  Average  error  (in  dB)  in  the  vertical-vertical  (vv)  analytical  solution  for  an 
untapered  circular  cylindrical  tree  trunk  as  functions  of  trunk  length,  trunk  radius,  and 
incidence  angle:  a)  6i  =  15°,  b)  6i  =  30°,  c)  6i  =  45°,  d)  6i  =  60°,  e)  6i  =  67°,  and  f)  6i  =  75° 
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a,  a, 

(e)  (0 


Fig.  3  Average  error  (in  dB)  in  the  horizontal-horizontal  (hh)  analytical  solution  for  an 
untapered  circular  cylindrical  tree  trunk  as  functions  of  trunk  length,  trunk  radius,  and 
incidence  angle:  a)  6i  =  15°,  b)  6i  =  30°,  c)  6i  =  45°,  d)  6i  =  60°,  e)  6i  =  67°,  and  f)  6i  =  75° 
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- 0.  =  15° 


Fig.  4  The  3-dB  error  lines  for  an  untapered  circular  cylindrical  tree  trunk  as  a  function 
of  incidence  angle:  vv  solution 


Fig.  5  The  3-dB  error  lines  for  an  untapered  circular  cylindrical  tree  trunk  as  a  function 
of  incidence  angle:  hh  solution 

Figures  2  and  3  show  that,  for  a  given  trunk  radius,  the  error  in  the  analytical 
solution  in  general  decreases  with  increasing  trunk  length  for  both  polarizations  and 
at  all  incidence  angles — albeit  the  error  does  not  decrease  in  a  monotonic  fashion, 
especially  for  small  length  values,  as  evident  by  the  non-smoothness  of  the  intensity 
maps.  The  error  profde  of  the  vv  solution  as  a  function  of  the  incidence  angle  is 
noticeably  different  from  that  of  the  hh  solution  owing  to  the  Brewster  angle  effects. 
The  range  of  validity  of  the  hh  solution  is  seen  to  improve  as  the  incidence  angle 
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increases,  but  there  is  an  incidence  angle  beyond  which  the  error  maps  remain 
relatively  unchanged;  here,  that  angle  seems  to  occur  at  around  =  45°.  In  contrast, 
the  vv  solution  tends  to  have  larger  errors  at  and  near  incidence  angles  at  which  the 
Brewster  angle  effects  of  the  trunk  and  the  ground  come  into  play.  In  this  case,  the 
Brewster  angles  of  the  trunk  and  the  ground  are  encountered  at  ~  15°  and  0i  ~ 
67°,  respectively,  and  it  is  seen  from  Figs.  2a  and  2e  that  the  vv  solution  accuracy 
is  degraded  at  those  angles. 

It  should  be  noted  that,  in  general,  as  the  single  ground-bounce  returns  become 
stronger  (or  more  dominant)  relative  to  the  other  scattering  effects,  the  accuracy  of 
the  analytical  solution  tends  to  improve.  This  can  be  explained  by  the  fact  that  the 
single  ground-bounce  responses  essentially  capture  the  scattering  behavior  of  the 
trunk  in  the  specular  direction  when  the  trunk  is  excited  either  directly  by  the 
incident  wave  or  indirectly  by  a  ground  reflected  wave. 

From  separate  free  space  scattering  simulations  of  the  finite-length  cylinder,  it  is 
seen  that  the  volumetric  current  integration  method  is  more  accurate  in  modeling 
the  specular  scattering  component  than  in  capturing  other  effects  such  as  scattering 
in  the  non-specular  directions  and  scattering  from  the  ends  of  the  structure. 
(Specifically,  the  accuracy  in  the  specular  direction  is  better  for  hh  than  for  w,  and 
improves  with  increasing  length  and  incidence  angle — with  the  exception  that  the 
vv  case  exhibits  an  error  peak  in  the  vicinity  of  the  Brewster  angle.)  Therefore,  in 
the  presence  of  the  ground,  if  strong  single  ground-bounce  returns  are  established, 
then  it  can  be  expected  that  the  analytical  solution  would  be  accurate.  For  the  vv 
response  at  6i  =  15°,  the  poor  accuracy  of  the  analytical  solution  is  due  to  the 
weakening  of  the  single  ground-bounce  returns,  the  specular  reflections  from  the 
trunk  not  being  faithfully  captured,  and  the  fact  that — especially  as  the  radius 
increases — the  backscattering  component  from  the  top  and  base  of  the  trunk 
becomes  stronger  than  or  comparable  to  the  single  ground-bounce  returns. 

Note  that,  at  least  for  vv,  the  analytical  formulation  better  characterizes  the 
backscattering  component  from  the  top  and  base  of  the  trunk  for  small  radius  values 
than  for  large  radius  values;  consequently,  at  6i  =  15°,  the  error  increases  with 
increasing  radius  for  a  given  length;  however,  at  other  incidence  angles,  on  average, 
the  error  tends  to  be  smaller  for  larger  radius  values  since  the  single  ground-bounce 
returns  are  the  more  dominant  responses  relative  to  other  effects.  As  it  can  be  seen 
by  juxtaposing  Figs.  2a  and  2e,  the  Brewster  angle  effect  of  the  ground  is  not  as 
detrimental  as  that  of  the  trunk;  at  0i  =  67°,  although  the  ground  reflection  is 
reduced,  the  finite-length  cylinder  approximation  still  better  captures  the  specular 
reflections  from  the  trunk  at  this  angle  than  at  0i  =  15°. 
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Overall,  the  error  behavior  of  the  hh  solution  exhibits  less  complexity  as  compared 
to  that  of  the  vv  case.  The  hh  backscattering  solution  has  less  than  3  dB  of  error  as 
long  as  (Fig.  5),  for  the  simulation  parameters  and  range  of  radius  and 

excitation  angle  values  considered  here.  The  stated  criterion  indicates  the  condition 
for  which  the  hh-polarized,  single  ground-bounce  returns  are  expected  to  be 
accurately  modeled.  It  is  also  observed  that  in  general  the  hh  analytical  solution 
may  not  always  accurately  account  for  the  scattering  from  the  top  and  base  of  the 
trunk,  even  for  small  radius  values. 

The  usefulness  of  the  error  analysis  conveyed  by  Figs.  2  through  5  notwithstanding, 
it  is  also  instructive  to  examine  the  scattering  characteristics  in  the  imaging  domain. 
The  imaging  response  at  pixel  location  f ,  for  a  particular  elevation  observation 
angle  6i,  is  obtained  by  averaging  the  backscattered  responses  over  an  aperture  in 
azimuth  as  follows: 


F™(F,6»,)  =  4;r 


1 


J^k^ki-r 


A/ 


,pp  =  w,hh 


(4) 

where  is  the  extent  of  the  observation  aperture,  is  the  number  of  observation 
angles,  )  is  the  window  in  azimuth,  and  Wf(f)  is  the  window  in  the  frequency 
domain.  In  this  work,  only  the  imaging  response  on  the  ground  plane  is  considered. 

Comparisons  of  the  FDTD  solution-based  and  analytical  solution-based  images  for 
a  few  representative  examples  are  shown  in  Figs.  6  and  7,  for  vv  and  hh, 
respectively.  These  examples  are  chosen  from  regions  where  the  average  error  of 
the  analytical  solution  in  the  frequency  domain  is  expected  to  be  less  than  3  dB. 

In  all  cases,  the  observation  aperture  is  centered  at  ^  =  0°,  with  ~  49°, 
=49 ,  A/=  [200  MHz,  500  MHz],  and  A^=  301;  and  the  base  of  the  trunk  is 

centered  at  the  origin. 
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Fig.  6  Imaging  response  (vv)  of  an  untapered  circular  cylindrical  tree  trunk, 
full  -wave  solution  (left)  vs.  analytical  solution  (right):  a-b)  ft  =  30°,  4=6.1,  =0.2; 

c-d)  di  =  60°,  =  3.6  ,  =  0.2 ;  e-f)  6»,  =  67°,  =  9.0  ,  =  0.3  ;  and  g-h)  ft=  75°, 

=  10.0,  =  0.5 
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Fig.  7  Imaging  response  (hh)  of  an  untapered  circular  cylindrical  tree  trunk, 
full-wave  solution  (left)  vs.  analytical  solution  (right):  a-b)  Qi  =  30°,  -  6.1,  -  0.2; 

c-d)  (9i=  60°,  lx,  "  ^-6  ,  =  0-2 ;  e-f)  (9,  =  67°,  Ix^  =  9.0  ,  =  0.3  ;  and  g-h)  (9,  =  75°, 

=  10.0,  ax  =  0.5 
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As  seen  from  the  results,  the  most  identifiable  features  in  the  images  come  from  the 
following:  1)  the  scattering  from  the  top  of  the  trunk,  2)  the  single  ground-bounce 
(dihedral)  returns,  and  3)  the  double  ground-bounce  response  of  the  top  of  the  trunk 
(that  is,  the  subsurface  mirror  term  of  the  scattering  from  the  top  of  the  trunk).  These 
are  referred  to  as  the  lower-order  effects,  with  the  single  ground-bounce  terms 
generating  the  dominant  effect  as  the  length  of  the  trunk  increases.  Note  that  the 
scattering  from  the  base  of  the  trunk  produces  imaging  responses  that  coincide  with 
the  single  ground-bounce  responses.  Of  course,  there  are  higher-order  (less 
prominent,  less  readily  identifiable)  effects  embedded  in  the  scattering  data:  these 
consist  of  higher-order  ground  bounces,  as  well  as  effects  produced  by,  for 
example,  internal  bounces  within  the  cylinder,  multiple  end-to-end  bounces, 
creeping  waves,  and  various  combinations  of  interactions  among  the  above  modes. 
It  can  be  shown  that  the  down-range  locations  of  the  3  lower-order  scattering 
features  are  +l^  /tanO. ,  x,  »  ,  and  ^3  «  -l^  Jtaae, ;  these  locations 

are  marked  in  the  images  by  dotted  vertical  lines  from  right  to  left. 

The  images  demonstrate  that  the  single  ground-bounce  returns  are  accurately 
modeled  by  the  analytical  solution;  the  scattering  from  the  top  of  the  trunk,  the 
double  ground-bounce  return,  and  higher-order  effects,  on  the  other  hand,  may  not 
be  captured  as  consistently,  or  as  accurately.  In  general,  considering  the  imaging 
response  as  a  whole,  the  match  between  the  numerical  and  the  analytical  solutions 
tends  to  improve  as  the  trunk  length  increases. 

3.2  Linearly  Tapered  Circular  Cylindrical  Trunk 

For  a  cone-shaped  tree  trunk  (Fig.  lb),  the  trunk-ground  dihedral  effect  is  reduced. 
In  other  words,  the  single  ground-bounce  returns  do  not  contain  the  direct  specular 
reflection  interactions  with  the  trunk  that  define  the  error  behavior  of  the  untapered 
case.  However,  when  the  cone  length  is  large  relative  to  the  base  radius,  the 
interactions  with  the  structure  are  such  that  a  scattering  component  in  a  direction 
close  to  specular  is  established  in  the  single  ground-bounce  returns. 

Error  maps  are  first  generated  as  before;  however,  for  space  considerations,  these 
maps  are  not  explicitly  included  here.  Instead,  only  their  3-dB  error  lines  are 
presented,  as  displayed  in  Figs.  8  and  9.  Comparing  Figs.  5  and  9  for  the  hh 
response,  it  is  seen  that  the  error  behaviors  for  the  untapered  and  linearly  tapered 
geometries  do  not  differ  significantly.  The  most  noticeable  difference  between  the 
2  cases  is  that,  for  a  fixed  length  value,  the  error  tends  to  increase  with  increasing 
base  radius  for  the  cone.  This  behavior — ^more  evident  for  small  length  values — is 
attributable  to  the  fact  that  as  the  base  radius  increases,  the  single  ground-bounce 
returns  contain  interactions  with  the  cone  that  are  farther  from  specular. 
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Fig.  8  The  3-dB  error  lines  for  a  linearly  tapered  circular  cylindrical  tree  trunk  as  a 
function  of  incidence  angle:  vv  solution 


Fig.  9  The  3-dB  error  lines  for  a  linearly  tapered  circular  cylindrical  tree  trunk  as  a 
function  of  incidence  angle:  hh  solution 

The  influence  of  the  Brewster  angle  effects  of  the  trunk  and  the  ground  on  the  w 
solution  error  is  more  complicated  than  that  observed  for  the  untapered  geometry, 
since  here  the  slant  angle  of  the  side  of  the  trunk  is  dependent  on  the  length-to-base 
radius  ratio.  Relative  to  the  Brewster  angle  effects,  the  single  ground-bounce 
returns’  non-specular  interactions  with  the  trunk  are  seen  to  be  the  more  important 
factor  in  determining  the  validity  of  the  analytical  solution.  As  a  result,  as  illustrated 
in  Fig.  8  (for  6i  =  30°,  45°,  60°,  and  75°),  the  region  of  validity  of  the  solution  for 
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the  cone  is  reduced  compared  to  that  for  the  untapered  cylinder;  the  error,  for  a 
fixed  length  value,  also  mostly  increases  with  increasing  radius  for  the  same  reason 
mentioned  above  for  the  hh  case.  Note  that,  relative  to  the  single  ground-bounce 
returns,  the  scattering  component  from  the  top  of  the  cone  is  weaker  than  the  one 
from  the  top  of  the  regular  cylinder,  as  expected.  For  the  vv  case,  this  leads  to  better 
accuracy  (larger  region  of  validity)  at  0i  =15°  and  67°  for  the  cone,  as  shown  in 
Fig.  8. 

Next,  images  are  generated  for  an  example  set  with  the  same  parameters  as  those 
considered  in  the  preceding  section  for  the  regular  cylinder.  The  results  are 
displayed  in  Figs.  10  and  11,  and  show  that  the  single  ground-bounce  components 
are  accurately  reproduced  by  the  analytical  solution.  Comparing  the  images  with 
those  for  the  regular  cylinder  case,  the  distinctive  feature  here  is  that  the  response 
from  the  top  of  the  structure  and  its  double  ground-bounce  return  have  been 
significantly  reduced  (relative  to  the  main  response);  this  is  reflected  in  both  the 
numerical  and  analytical  solutions.  (Note  that  it  is  not  appropriate  to  make  a  direct 
comparison  of  the  amplitudes  of  the  image  responses  between  the  regular  cylinder 
and  cone  cases  here,  since  for  the  same  length  and  radius  values,  the  2  structures 
have  different  volumes.) 
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Fig.  10  Imaging  response  (vv)  of  a  linearly  tapered  circular  cylindrical  tree  trunk, 
full-wave  solution  (left)  vs.  analytical  solution  (right):  a-b)  ft  =  30°,  =6.1,  =  0.2; 

c-d)  ft—  60°,  =  3.6  ,  =  0.2  ;  e-f)  ft—  67°,  =  9.0  ,  =  0.3  ;  and  g-h):  ft—  75°, 

L.  =  10.0,  a.  -  0.5 

4  ’  4 
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Fig.  11  Imaging  response  (hh)  of  a  linearly  tapered  circular  cylindrical  tree  trunk, 
full-wave  solution  (left)  vs.  analytical  solution  (right):  a-b)  Qi  =  30°,  -  6.1,  -  0.2; 

c-d)  (9i=  60°,  =  3.6  ,  =  0.2 ;  e-f)  (9,  =  67°,  =  9.0  ,  =  0.3  ;  and  g-h)  (9,  =  75°, 

=  10.0,  =  0.5 
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3.3  Nonlinearly  Tapered  Circular  Cylindrical  Trunk 


The  trank  geometry  examined  next  is  akin  to  a  cone  structure  with  an  exponential 
expansion  near  its  base.  The  tapering  profile  shown  in  Fig.  Ic  is  generated  by 
following  the  procedure  outlined  by  Weber  and  Penn''  using  the  parameters 
nTaper  =  1  and  Flare  =  0.6.  As  mentioned,  such  a  profile  with  a  slightly  flared  base 
is  expected  to  provide  a  closer  emulation  of  trank  geometries  encountered  in  nature. 

Only  the  3-dB  error  lines  for  the  analytical  solution  are  included,  as  shown  in 
Figs.  12  and  13.  The  plots  are  similar  to  those  for  the  cone  structure.  It  is  seen  that 
the  presence  of  the  flare  does  not  significantly  change  the  error  behavior.  The 
observations  presented  previously  for  the  cone  are  equally  applicable  here.  (Note 
that,  as  shown  by  Weber  and  Penn,"  by  changing  the  nTaper  and  Flare  parameters, 
a  variety  of  trank  profiles  can  be  generated;  the  values  used  here  are  expected  to  be 
adequate  for  the  types  of  trees  treated  in  this  study.) 

Images  for  tranks  with  the  above  tapering  profile  are  illustrated  in  Figs.  14  and 
15 — for  the  same  example  set  considered  in  the  previous  2  sections.  Good 
agreement  between  the  numerical  and  analytical  solutions  is  seen — though,  once 
again,  the  single  ground-bounce  returns  are  more  accurately  reproduced  than  other 
less  prominent  scattering  effects.  Overall,  the  results  demonstrate  that  the  discrete 
cylinder  model  can  be  an  acceptable  approach  for  analyzing  a  trank  with  non¬ 
constant  radius  along  its  length. 


Fig.  12  The  3-dB  error  lines  for  a  nonlinearly  tapered  circular  cylindrical  tree  trunk  as  a 
function  of  incidence  angle:  vv  solution 
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Fig.  13  The  3-dB  error  lines  for  a  nonlinearly  tapered  circular  cylindrical  tree  trunk  as  a 
function  of  incidence  angle:  hh  solution 
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Fig.  14  Imaging  response  (vv)  of  a  nonlinearly  tapered  circular  cylindrical  tree  trunk, 
full-wave  solution  (left)  vs.  analytical  solution  (right):  a-b)  Qi  =  30°,  -  6.1,  -  0.2; 

c-d)  (9i=  60°,  =  3.6  ,  =  0.2 ;  e-f)  <9,  =  67°,  "  9.0  ,  =  0.3  ;  and  g-h)  (9,  =  75°, 

=  10.0,  =  0.5 
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Fig.  15  Imaging  response  (hh)  of  a  nonlinearly  tapered  circular  cylindrical  tree  trunk, 
full-wave  solution  (left)  vs.  analytical  solution  (right):  a-b)  Qi  =  30°,  -  6.1,  -  0.2; 

c-d)  di  =  60°,  =  3.6  ,  =  0.2 ;  e-f)  6»,  =  67°,  =  9.0  ,  =  0.3  ;  and  g-h)  ft=  75°, 

=  10.0,  =  0.5 
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4.  Conclusions 


The  accuracy  of  the  volumetric  current  integration  method  for  characterizing  the 
backscattering  response  of  a  tree  trunk  located  above  a  dielectric  ground  is 
investigated.  The  region  of  validity  of  the  analytical  solution — as  functions  of  trunk 
length,  trunk  radius,  incidence  angle,  polarization,  and  trunk  tapering — is 
determined  by  comparing  the  frequency-domain  results  to  those  from  a  full-wave 
method.  It  is  noted  that  the  w  solution  is,  in  general,  valid  over  a  smaller  range  of 
length  and  radius  values  than  the  hh  solution,  as  a  result  of  the  Brewster  angle 
effects  of  the  trunk  and  the  ground.  The  range  of  validity  is  also  negatively 
impacted — albeit  less  so  for  hh  than  for  vv — by  the  presence  of  tapering,  a  feature 
that  results  in  single  ground  bounces  defined  by  non-specular  interactions  with  the 
trunk.  It  is  observed  that  for  both  the  untapered  and  tapered  profiles,  the  range  of 
validity  of  the  2  co-polarized  responses  becomes  comparable  as  the  incidence  angle 
approaches  grazing.  To  demonstrate  the  practicality  of  the  formulations  and 
observations,  monostatic  images  of  the  tree  trunks  are  also  constructed. 
Comparisons  with  numerical  solution-based  images  indicate  that  while  the 
analytical  solution  can  correctly  model  the  single  ground-bounce  returns,  it  may 
not  consistently  capture  other  less  prominent  effects  such  as,  for  example,  scattering 
from  the  top  of  the  structure  and  higher-order  interactions. 
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Appendix.  Scattering  Solution  for  the  Finite-Length,  Dielectric 

Circular  Cylinder  in  Free  Space 
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Formulations  from  Part  I  of  this  study^  are  included  here  for  convenience.  The  free 
space  scattering  matrix  elements  of  a  finite-length,  dielectric  circular  cylinder  with 
length  L,  radius  a,  and  complex  relative  dielectric  constant  Sr,t  are  derived  as 


{kA)  =  Ps 


COS  0.  +  cos  0^  ) 
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The  above  assumes  k^.  =k  sine,  k„  =k  sinO  ,  k  Js^^-cos^  a  ,  and  that 

yC^jZ  O  I  P  O  S  H  L/  ^  r  I 

/V  yv 

the  incident  and  scattered  wave  directions  ( k. ,  )  and  their  associated 

polarization  vectors  {h  ,  v)  are  defined  by 


=  sin  9.  cos  (p.x  +  sin  9.  sin  </).y  -  cos  9^1 . 
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4  =  sin  9^  cos  +  sin  9^  sin  +  cos  9J . 


zxk, 
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(A- 10) 
(A-11) 
(A-12) 

(A- 13) 


in  which  the  subscripts  i  and  s  identify  quantities  related  to  the  incident  and 
scattered  waves,  respectively;  0i  is  the  incidence  angle  in  elevation  (measured  from 
-  z  );  6s  is  the  scattering  angle  in  elevation  (measured  from  +  z  );  and  ^,s  is  the 
incidence/scattering  angle  in  azimuth  (measured  from  +  x). 
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